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Temperature-dependent resistivity is studied in single crystals ol iron-arsenide superconductor 
Nai-iFei-^Coa; As for electrical current directions along, p a (T), and transverse, p c (T), to the Fe-As 
layers. Doping with Co increases stability of this compound to reaction with the environment and 
suppresses numerous features in both p a {T) and p c (T) compared to the stoichiometric NaFeAs. 
Evolution of p a {T) with x follows a universal trend observed in other pnictide superconductors, 
exhibiting a T-linear temperature dependence close to the optimal doping and development of T 2 
dependence upon further doping. pc(T) in parent compound shows a non - monotonic behavior 
with a crossover from non-metallic resistivity increase on cooling from room temperature down to ~ 
80 K to a metallic decrease below this temperature. Both p a (T) and p c (T) show several correlated 
crossover - like features at T > 80 K. Despite a general trend towards more metallic behavior of 
inter - plane resistivity in Co-doped samples, the temperature of the crossover from insulating to 
metallic behavior (80 K) does not change much with doping. 

PACS numbers: 74.70.Dd,72.15.-v,74.25.Jb 



I. INTRODUCTION 

Structurally, iron based superconductors are layered 
materials, in which FeAs (or iron chalcogenide) layer is 
the main building block for a variety of compoundJ^l. 
Since the dominant contribution to the density of states 
at the Fermi level comes from the iron 3d orbitals, one 
can expect a significant electronic anisotropy of the com- 
pounds revealed in the in-plane and out-of-plane trans- 
port. Contrary to this expectation, the most studied 
families of iron arsenides, those based on BaFe2As2, have 
rather low anisotropy ratio 7 p = p c /p a ^4 at T^. In tran- 
sition metal-doped Ba(Fei_ a; M x ) 2 As 2 (M=Co, Ni, Rh, 
Pd, BaT122 in the following), p c (T) also shows a very 
different temperature dependence compared with p a (T), 
revealing a broad crossover from non-metallic to metal- 
lic temperature dependence assigned in our systematic 
doping studies to the formation of a pseudogajPtil. 

Another interesting feature of iron arsenides that dis- 
tinguishes them from the copper oxide based (cuprate) 
superconductora^ is a strong variation of the functional 
form of temperature-dependent resistivity for various 
types and levels of dopings. The general trend of p a (T) 
evolution is the presence of a T— linear region imme- 
diately above T c for optimally doped compositions^^. 
At higher temperatures this T-linear behavior, for ex- 
ample in hole-doped Bai_ a: K x Fe2As2 and in sel f-doped 
Nai-^FeAs, is terminated by the pseudogap F 2 * 13 l 

Systematic studies of the temperature-dependent elec- 
trical resistivity are very important for the general under- 
standing of superconductivity in this family of materials. 
Scattering in the normal state in the vicinity of the mag- 
netic quantum critical point leads to the characteristic 
T— linear temperature dependence of p a (T), which then 
evolves towards Fermi-liquid T 2 — behavior with doping 
(see Ref. [T3] for a review). Deviations from this general 



behavior provide an insight into electronic and magnetic 
correlations^, in particular, into the mechanism of ne- 
matic state formatiorP^lEl] 

In this article we report the systematics of doping- 
evolution of in-plane and inter-plane resistivity of 
electron-doped Nai-sFei-zCo^As. This compound 
shows a "d ome - lik e" phase diagram which is very similar 
to BaCol2 2 12 l 19 l 2t H As such, this study brings additional 
insight into the scattering and correlation phenomena of 
the iron - based superconductors. 



II. EXPERIMENTAL 

AC magnetic characterization of the samples was 
performed with a tunnel-diode resonator, (TDR)P^U. 
Briefly, TDR is a self-oscillating LC tank circuit pow- 
ered by a properly biased tunnel diode. The sample is 
mounted with Apiezon N-grease on a sapphire rod and 
is inserted in the inductor (coil). The sample temper- 
ature is controlled independent of the resonant circuit, 
which is actively stabilized at the constant temperature. 
The measured frequency shift is proportional to the dif- 
ferential magnetic susceptibility of the sample^. In this 
work, for quick mounting and measurement protocols we 
used a simplified version of the TDR susceptometer (a 
"dipper" ) , which is inserted directly into a transport He 
dewar and gives very quick turn-around measurement 
time of typically 30 minutes per sample. The trade-off 
of this quick measurement protocol is reduced stability 
and higher temperature-dependent background as com- 
pared to our high-stability 3 He and dilution refrigerator 
versions of the TDR susceptometer. Nevertheless, the 
"dipper" is perfectly suitable to study magnetic signa- 
ture of the superconducting transition. 

Single crystals of Na(Fei_ x Co :E )As with x = 0, 0.025, 
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0.05, 0.08, and 0.10 were synthesized by sealing a mix- 
ture of Na, Fe, As, and Co together in Ta tubes and heat- 
ing it to 950 °C, followed by 5 °C/hour cooling down to 
900 °CP3. The x in samples was defined as the nomi- 
nal ratio, which give s som e variation with electron-probe 
microanalysis value^2122l The samples were stored and 
transported in sealed containers filled with inert gas. 

Sample preparation was done quickly in air within 
about 5 minutes to minimize uncontrolled environmen- 
tal exposure which can induce an increase of TtPHU. We 
started sample preparation by cleaving slabs from the 
inner part of the crystals with typical a thickness of 50 
to 100 pm. The slabs had shiny cleavage surfaces and 
were further cut into smaller pieces for TDR (typically 
0.5 x 0.5mm 2 ) and resistivity measurements. Cleaved in- 
ternal parts of single crystals did not show any visible 
reaction with air and turned out to be relatively stable, 
contrary to crystals with the residue of NaAs flux, which 
aggressively reacts with air and moisture. After prepara- 
tion, samples were promptly measured and immediately 
stored after measurements in inert and dry environment. 
After each dipper run samples were washed with toluene 
to remove remaining N-grease in order to control the air 
exposure. 

Samples for in-plane resistivity measurements had typ- 
ical dimensions of (l-2)x0.5x(0. 02-0.1) mm 3 . All sam- 
ple dimensions were measured with an optical micro- 
scope with an accuracy of about 10%. Sample resis- 
tivity at room temperature, p(300-ftT), was in the range 
400 to 500 pilcm for all compositions studied. This 
value is obtained on a bigger array of samples than in 
our previous study-^ and is somewhat higher. It is also 
somewhat hig her than values found in electron-^21 anc j 
hole-doped^ Bal22 compounds, typically 300 pVLcm or 
less. We do not have sufficient array of data to obtain 
lower error bars needed to resolve the doping evolution of 
p(300K), if any exists. Contacts for four-probe resistivity 
measurements were made by soldering 50 pm silver wires 
with ultra-pure Sn solder, as described in Ref. [3H] . Resis- 
tivity measurements were performed in QuantumDesign 
PPMS, providing magnetic fields up to 9 T. For mea- 
surements of the upper critical field, H C 2, samples were 
glued to the side of a plastic block with ab plane of the 
sample oriented to be either parallel or perpendicular to 
the direction of magnetic field (with an accuracy of about 
1°). 

Inter-plane resistivity measurements were done using 
the two-probe technique, relying on very low contact 
resistance of soldered contacts, typically in the 10 pQ, 
range. The top and bottom surfaces of the a6-plane 
(typically 0.5 x 0.5mm 2 area) of the samples were cov- 
ered with Sn solder forming a capacitor-like structure. A 
four-probe scheme was used to measure a sample with 
contacts, giving a sum of series connected sample, R s , 
and contact resistance, R c resistances. Since R s 3> R c , 
contact resistance represents a minor portion, on the or- 
der of 1-5% on the total resistance. This can be directly 
seen for our samples for temperatures below the super- 



conducting T e , wher e R s =0 and the measured resistance 
represents jjH 4 l 28 l 29 l. Further details of the measurement 
procedure can be found in Refs. 0HS1 

The drawback of the measurement of samples with 
c <C a is that any structural and chemical inhomogeneity 
along the c— axis, a very common problem in soft and mi- 
caceous samples of iron arsenide superconductors ^ 30 * 32 !, 
not only increases sample resistance, but admixes in- 
plane component due to the redistribution of the cur- 
rent. One way to ascertain correctness of the p c mea- 
surements, is to rely on measurements with the lowest 
resistivity values. Typically the best results were ob- 
tained on the thinnest slabs. To get reliable results we 
performed measurements of p c on at least 5 samples of 
each batch. In all cases we obtained qualitatively similar 
temperature dependences of the normalized electrical re- 
sistivity, p c (T) / p c (300K). The resistivity value at room 
temperature, p c (300K), however, showed a notable scat- 
ting and was typically in the range 2000 to 3000 pfl cm 
at room temperature. 

We have shown previously that reaction with air 
strongly affects the value of p a (300-fQ due to the develop- 
ment of cracks^. Cracks grossly effect the internal sam- 
ple connectivity and, hence, homogeneous current dis- 
tribution, thus, making inter-plane resistivity measure- 
ments of environmentally exposed samples impossible. 

III. RESULTS 
A. Environmental stability 
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FIG. 1. (Color online) Evolution of the frequency shift 
signal in the dipper TDR experiment on increasing time 
of air-exposure treatment in crystals of slightly overdoped 
Nai-aFei-^Co^As x=0.08. Similar to parent compouncpl 
T c of the sample increases with exposure time to a maximum 
and then decreases. Similar effects are observed for other 
doping levels, x. 

The T c of the parent Nai_^F eAs increases signifi- 
cantly upon exposure to aii ! 12 * 33 * , water^l and Apiezon 
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FIG. 2. (Color online) Doping evolution of the frequency shift 
signal in the dipper TDR experiment during fixed time air- 
exposure treatment of NaFei-^Co^ As for 24 hours. As can be 
seen, T c increase from 12 K to 22 K, characteristic of parent 
NaFeAs samples, is strongly suppressed with Co doping. 
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signal during fixed time, one day air exposure, is summa- 
rized in Fig. [2j Variation of fresh sample T c and highest 
achieved T c during air exposure and fixed - time of one 
day exposure as a function of x are summarized in Fig. [3] 

One day exposure of a sample to air does not lead to 
a visual appearance of reaction products. Thus, at least 
at this initial stage, there is no reason to assume trans- 
formation of NaFeAs into NaFe2As2, the final product 
of reaction with watei^, formed after one a month ex- 
posure, which most likely shows up in Fig. [I] as a new 
shoulder in the temperature dependent frequency shift 
at about 12 K for samples exposed for about two weeks. 

In parent Nai^FeAs the environmental reaction is 
caused by the variation of Na content in the samples, 
6, due to oxidative deintercalatiorP^. It is natural to 
expect similar effect in the Co - doped NaFei-^Co^As. 
However, the puzzling observation is that T c increases 
for both environmental reaction in pure Nai_,$FeAs (pre- 
sumably hole - like doping) and electron Co - doping. As 
such, it is not clear if carrier type and density change is 
the main effect involved. Wc note that detailed study 
of the effect of Li deficiency in a closely related LiFeAs 
superconductor found suppression of Tc, but virtually no 
change in the normal - state propertied^. If Na deficiency 
leads to the formation of Na vacancies, this should lead 
to hole doping and, thus, move the dome on the doping 
phase diagram in an opposite way to electron Co-doping. 
Further studies are required to understand what type of 
doping is induced by the loss of Na and what types of 
defects are formed. 

At a first glance, the different rate of T c variation in 
Fig. [3] can be attributed to the different sensitivity of 
T c (x) in different parts of the phase diagram - being 
smallest at the flat optimal doping region. Indeed, the 
slope of T c (x) changes from very high for the parent com- 
pound to negligible in samples close to optimal doping. 
However, we find a rise in T c after environmental reaction 
even in over-doped samples. This contradicts the simple 
relation of the rate of T c change to be determined by a 
position on the phase diagram. Our observations rather 
support previous observation 22 that Co doping increases 
the stability of the samples. 



FIG. 3. (Color online) Doping evolution of the superconduct- 
ing Tc in the dipper TDR experiment during air-exposure 
treatment of NaFei-^CoaAs. Blue solid dots show T c of fresh 
samples, red down-triangles T c of the samples exposed to air 
for one day. Green up-triangles show maximum onset T c ob- 
tained in these experiments. 



N-grease^. Here, we study how the sensitivity of T c to 
exposure changes with Co- doping. In Fig. [T] we show 
the evolution of the TDR signal in slightly overdoped, 
£=0.08 (fresh sample T c =16 K), samples on exposure 
to air. Similar to the parent compound^, for all com- 
positions irrespective of their x, the T c of the samples 
increases initially upon air exposure and then decreases 
with prolonged exposure. The doping-variation of TDR 



B. Resistivity measurements 

The temperature-dependent resistivity of "fresh" crys- 
tals of NaFei-zCo^As is shown in Fig. [I using a nor- 
malized resistivity scale, p/ p(300K). The shape of p a {T) 
in the parent compound is relatively complex, with fea- 
tures due to split structural (at temper ature T s =55 K) 
and magnetic (at T m =45 K) t ransit ion d 19 * 36 ! and slope 
changes at higher temperature d 12 * 3 ^ . With doping the 
dependence transforms to very close to T-linear p a (T) in 
samples with £=0.08 and close to T 2 on further x in- 
crease. The changes of slope at T\ ^300 K (increase of 
slope on cooling), T 2 ^160 K (decrease of slope on cool- 
ing) , and T3 ~80 K (increase of slope on cooling) are ob- 
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FIG. 4. (Color online) Doping-evolution of the temperature 
dependence of the in-plane resistivity, p a /pa(300K), for sam- 
ples of NaFei-ajCo^As in fresh state after initial sample han- 
dling and contact making. For reference we show data from 
Ref. [SUfor stoichiometric LiFeAs, representative of the over- 
doped regime. 
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FIG. 5. (Color online) Doping-evolution of the temperature 
dependence of the inter-plane resistivity, p c / p c (300K), for 
samples of NaFei-^Co^As in fresh state after initial sample 
handling and contact making. For reference we data from 
Ref. [34|fbr stoichiometric LiFeAs, representative of the over- 
doped regime. 



served in doped samples, similar to the parent compound, 
and are relatively insensitive to doping. The feature at 
T\ is observed in samples with all x studied. It is similar, 
though less pronounced, to a slope change at about the 
same temperature in p a (T) of stoichiometric LiFeAs, see 
Fig. [4j The features at T 2 and T 3 are observed in p a (T) 
of the samples with x < 0.05. 

The results of this study of p a {T) are in reaso nable 
agreement with previous studies on sing le crystalPM] 
with the difference of x coming from using nominal values 
during sample preparation. The doping-transformation 
of the temperature dependent resistivity for T right 



above T c follows general expectations for a quantum crit- 
ical scenario^, with T-linear range confined from high 
temperature side by slope change on approaching T% . By 
comparison with position of maximum in the tempera- 
ture dependent inter-plane resistivity p c (T), we assigned 
similar slope change feature in p a (T) of Bai_; C K. r Fe2As2 
to formation of pseudogap^. The slope changes upon 
cooling through T 2 and T3 in samples with x <0.05 do 
not have a direct analogy with Bal22 compounds. These 
features are observed even in samples in which long-range 
magnetic order and orthorhombic structural distortion 
are suppressed. Studies of resistivity anisotropy on de- 
twinned single crystals of parent NaFeAs^ suggest that 
the feature at T3 has a similar nature to nematic cor- 
relations , whi ch is particularly strong in electron doped 
BaCol223Sini. 

For understanding the resistivity of NaFeAs based 
compounds, it is important to get insight into the tem- 
perature dependence of the inter-plane resistivity com- 
ponent. In Fig. [5] we show the doping evolution of p c (T) 
in NaFei-^Coz As. The inter-plane resistivity of parent 
NaFeAs increases during cooling down to a maximum at 
~70 K, which is close to T3 as determined from p a (T). 
The resistivity rapidly decreases below this maximum, 
with a notable rate increase below T m . Note that con- 
trary to p a (T), the inter-plane resistivity does not show 
an increase below T s , suggesting that the carriers affected 
by the formation of a gap do not contribute much to the 
inter-plane transport. Interestingly, despite the strong 
difference between p a (T) and p c (T), the high tempera- 
ture features are observed in both of them. 



C. Anisotropy of the upper critical field 
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FIG. 6. (Color online) Evolution of the temperature- 
dependent resistivity in the vicinity of the superconducting 
transition in magnetic fields applied perpendicular to (H || c, 
left panel) and parallel to (H \\ ab, right panel) the conducting 
Fe-As plane of the sample of optimally doped NaFei-^CoajAs, 
x=0.025. The resistive transition temperature, used to plot 
H — T phase diagram shown in Fig.[7]below, was defined using 
midpoint criterion. 

The anisotropy of the electrical resistivity at T c , 7 P = 
^ c ^ c j , is linked with the anisotropy of the upper critical 
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FIG. 7. (Color online) Temperature- magnetic field phase 
diagram of optimally doped NaFei-^Co^As, £=0.025 for ori- 
entation of magnetic field perpendicular and parallel to the 
conducting Fe-As plane of the crystal. 

field, 7# = ^gr^wr > with 7p = 1h- Because determi- 
nation of the absolute values in resistivity measurements 
always includes uncertainty of the geometric factor and 
is affected by the cracks, the H C 2 anisotropy measure- 
ments provide an alternative way to evaluate resistivity 
anisotropy^. In Fig. [6] we zoom the superconducting tran- 
sition in in-plane resistivity measurements p a {T) for a 
sample with doping level close to optimal, £=0.025. The 
same sample was remounted on a plastic cube with the 
magnetic field was in the H || c (left panel) and H || ab 
(right panel) configurations. 

We used the resistive transition midpoint to determine 
H C 2(T) anisotropy as shown in Fig. [7] Close to T c the 
anisotropy 7#=2.25 ±0.1 for the sample with £=0.025. 
Similar value with 7^=2.35 ±0.1 was obtained in sam- 
ple with £=0.08. These measurements suggest a resis- 
tivity anisotropy of about 5 at T c . Considering that 
7 P (T C ) ~ 27 p (300 K) (see Figs. [2] and |) we expect a 
negligible anisotropy of 2 to 3 at room temperature. The 
direct resistivity measurements, with p a (300i4T)=400 to 
500 pttcm and p c (300X)=2000 to 3000 pflcm, suggest 
an anisotropy of 4 to 8. The origin of this factor of about 
two discrepancy remains unclear at the moment. 



IV. DISCUSSION 

A. Slope-change features in the 
temperature-dependent resistivity 

As can be seen from direct comparison of in-plane and 
inter-plane resistivity in the parent and slightly doped 
£=0.025 compositions, Fig. [8j features in p a (T) find 
counterparts in p c {T). For example, a slope decrease 
in generally metallic p a (T) below T2 ^160 K is seen as 
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FIG. 8. (Color online) Comparison of the temperature- de- 
pendent in-plane and inter-plane resistivity of parent NaFeAs 
(left panel) and a sample with a;=0.025 (right panel). Resis- 
tivity data are plotted vs. normalized scale, p(T)/p(300A"). 
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FIG. 9. (Color online) Comparison of the temperature- de- 
pendent in-plane resistivity in samples of NaFeAs doped to 
the highest T c by environmental reaction with Apiezon N- 
grease and with Co-doping, £=0.05. Resistivity data are plot- 
ted vs. normalized scale, p(T) / p(300K). Lines show linear 
extrapolation of the curves to T — >0, revealing difference in 
residual resistivity of two concentrations of samples. 



slope increase in generally activated (increasing on cool- 
ing) p c (T). This similarity found in two very different 
dependences suggests that the activation of carriers over 
a partial gap, rather than change of scattering, is respon- 
sible for the feature. Partial (nematic) order^, which 
happens above the structural transition at T3, changes 
p c (T) from insulating to metallic, while the magnetic or- 
der below T m causes a dramatic decrease of resistivity in 
both directions of charge flow. The decrease is especially 
strong in the parent compound in which the residual re- 
sistivity ratio (RRR), p(300K )/p(T c ), is a factor of two 
higher than in £=0.025. These observations suggests that 
magnetic scattering plays an important role in resistivity 
at T > T3 , and that when inelastic scattering is dominant 
(as in the parent compound) taming down of magnetic 
fluctuations reveals intrinsically very low residual resis- 
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tivity. 



B. Residual resistivity 

Observation of a much higher RRR in non-doped ma- 
terials agrees with studies in Bal22 compounds, though 
in the latter, the direct comparison is not so simple. In 
the case of NaFeAs based materials we can compare RRR 
of the samples, brought to optimal doping using two dif- 
ferent doping, electron with Co substitution of Fe and 
environmental, on interaction with the environment. In 
Fig.[9]wc compare T-dependent resistivity in two concen- 
trations of samples, extrapolating curves linearly from T c 
to T = 0. The RRR ratio decreases from more than 20 in 
environmentally doped samples to about 4 in Co-doped 
samples. Taking that resistivity at room temperature 
does not change from about 400 /xficm , this suggests 
that the residual resistivity induced by the x=0.025 sub- 
stitution of Fe atoms with Co is on the order of 100 /iSlcm, 
comparable to BaCol22 e . This is almost a factor of five 
higher than the value found in very disordered samples, 
doped with environmental reaction, which extrapolates 
to po «20 pflcm. 
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